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The conformal gravity is one of the most important models of quantum gravity with higher 
derivatives. We investigate the role of the Gauss-Bonnet term in this theory. The coincidence limit 
of the second coefficient of the Schwinger-DeWitt expansion is evaluated in an arbitrary dimension 
n. In the limit n = 4 the Gauss-Bonnet term is topological and its contribution cancels. This 
cancellation provides an efficient test for the correctness of calculation and, simultaneously, clarifies 
the long-standing general problem concerning the role of the topological term in quantum gravity. 
For n / 4 the Gauss-Bonnet term becomes dynamical in the classical theory and relevant at the 
quantum level. In particular, the renormalization group equations in dimension n — 4 — e manifest 
new fixed points due to quantum effects of this term. 
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I. INTRODUCTION 

At both classical and quantum levels, local conformal 
symmetry plays a special role in theories of gravity and 
their applications [1-3]. One of the most interesting is- 
sues is the violation of this symmetry at the quantum 
level. For the quantum theory of matter fields in curved 
space-time the violation of conformal symmetry is re- 
lated to the well-known trace anomaly (see, e.g. [1, 3] 
for the review). The important feature of the conformal 
anomaly is its universality for the matter (scalar, spinor 
and vector) fields which contribute with the same sign 
to the two of three terms of the anomaly. The opposite 
sign of the contributions takes place for the unphysical 
higher derivative scalars and fermions [4, 5]. In princi- 
ple, one can choose the number of these higher derivative 
fields in such a way that they cancel the contributions of 
the matter fields. In this case the conformal symmetry 
holds at the one loop level. The cancellation of anomaly 
can not be achieved in the known versions of conformal 
supergravity [6], and therefore the relation between the 
cancellation of conformal anomaly and what is supposed 
to be the fundamental theory (e.g. supergravity which 
may be a low-energy limit of the (super)string/M - the- 
ory) remains unclear within the semiclassical approach. 

The violation of the conformal symmetry in quantum 
gravity is much less studied. One of the simplest theories 
of gravity which possesses local conformal symmetry is 



based on the Weyl action j d A 



where 
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is the square of the Weyl tensor in n = 4 dimensions. 
In order to provide renormalizability, one has to include 
topological and surface terms. In this way we arrive at 
the action 
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J d 4 xy/^j {^C 2 + 7]E + T OR } . 



Here 



E = R 



R» 



AR a pR 



a/3 



R' 



(2) 



(3) 



is the integrand of the Gauss-Bonnet topological term. 
The action (2) is conformal invariant, for it satisfies the 
conformal Noether identity 
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By dimensional reasons one can introduce into action 
(2) an extra term ■ J y/—g R 2 . However, this expres- 
sion possesses only global and not local conformal sym- 
metry and hence it will not be included into the action. 
In order to complete the story, let us notice that a fi- 
nite 9 ■ J \J—g R 2 term may be generated as a quantum 
anomaly-induced correction, e.g. due to the renormaliza- 
tion of the J yJ—gE\R term in (2) . The anomalous viola- 
tion of local conformal symmetry in the finite part of the 
one-loop effective action may produce the non-conformal 
divergences beyond the one-loop level. This effect has 
been investigated in [7] for the conformal scalar field and 
there are no reasons to expect that the situation for the 
conformal quantum gravity will be different. 

The main purpose of the present paper is the one loop 
renormalization and renormalization group in the con- 
formal quantum gravity. The renormalization structure 
depends on the form of divergences and corresponding 
counterterms. According to the standard expectations, 
despite the anomaly results from the one-loop renormal- 
ization, one-loop divergences in conformal quantum grav- 
ity must be conformally invariant. This property holds in 
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all known examples of conformal matter fields, and one 
can expect that the same should be true for the Weyl 
quantum gravity based on action (2). The most natural 
result would be to meet the renormalization of the coeffi- 
cients r], A, t but not the J \J~^g R 2 -type counterterm. 

From the first sight there is no much difference whether 
the / ^J~ r g~R 2 term shows up already in the one-loop di- 
vergences or only at higher loops. However, this may 
be relevant for some applications of quantum gravity. If 
non-conformal divergences show up only at higher loops, 
the conformal symmetry may be considered as a good 
approximation. For example, the one-loop renormaliz- 
ability of conformal gravity provides the possibility of 
the successful realization of the anomaly-induced infla- 
tion scheme (see the discussion in [8]) in the presence of 
quantum gravity. At the same time, if the non-conformal 
divergence emerges at the one-loop order, the conformal 
symmetry can not be considered a reasonable approxi- 
mation, because the running of the coefficient 9 will be 
much stronger in this case. In this situation a quantized 
conformal matter on curved classical background also can 
not be considered as an approximation to the full theory 
involving quantum gravity (see, e.g. discussion in [1, 9]). 
Finally, we need to be sure whether the non-conformal 
divergence is present at the one loop level in the Weyl 
quantum gravity. 

The first explicit derivation of the one- loop divergences 
in Weyl quantum gravity has been performed by Frad- 
kin and Tseytlin [10] in the framework of background 
field method, properly modified for the higher derivative 
theories [31]. The J y/—gR 2 -type divergence has been 
encountered and the conformal invariance of the coun- 
terterms has been achieved through the use of the special 
procedure of conformal rcgularization. This regulariza- 
tion is nothing but the specific rcparametrization of back- 
ground metric invented earlier in [13] (see also [14] and 
[15]). According to this procedure metric g^ v has to be 
replaced by conformal metric <7 M „ = g^ v P 2 [g M1/ ], where 
scalar metric-dependent quantity P[<7 M „] is defined as a 
solution of the equation 

OP=-RP. (5) 
6 v ' 

When performing a local conformal transformation of 
the original metric 

-> d'ixv = 9\u> e 2a{x) , the quantity 
P[g^ v ] transforms as P -> P' = Pe^ (l) , 

such that the metric g^ v remains invariant. Another 
important property of the metric g^ v is that the cor- 
responding scalar curvature is zero R = Rig^u) = 0. 
Therefore, after the original metric g^ v is replaced by 
g^v, the divergent j y 1 —gR 2 counterterm disappears 
and the expected invariant form of divergences gets re- 
stored. The procedure of conformal "regularization" has 
been generalized for the conformal quantum gravity cou- 
pled to conformal quantum matter fields in [15]. 



Is it correct to consider the replacement g^ — > g^ as 
a kind of conformal regularization for the divergent part 
of the effective action of quantum gravity? It is easy to 
see that this procedure eliminates also the anomaly in 
the finite part of the effective action [32]. Therefore, this 
choice of background metric does not fit with numerous 
applications of conformal anomaly which we know. Fur- 
thermore, despite the choice of g^u as a background met- 
ric is mathematically consistent, it is not very appealing 
because, in particular, it eliminates the Einstcin-Hilbcrt 
action. Consequently, the theory based on this metric 
may not have a consistent non-rclativistic limit. In gen- 
eral, the whole procedure looks as an artificial addition to 
the background field method. If we really want to learn 
the role of the conformal symmetry in quantum gravity, 
it is important to know whether the appearance of the 
J ^J—gR 2 counterterm is a calculational error or it is 
caused by inconsistency of the background field method 
applied to Weyl quantum gravity. The last option has 
been partially explored in [16], where the possible con- 
flict between diffeomorphism and conformal gauge fix- 
ing conditions has been discussed. It turned out that 
counterterm j y '—gR 2 is gauge fixing independent, ex- 
actly as the renormalization of the terms J \/—gC 2 and 
J ^J~^gE. On the opposite, counterterm J ^y—gOR is 
not protected from the gauge fixing dependence and can 
be modified or even eliminated by the appropriate choice 
of the parameters of gauge fixing. 

The second explicit derivation of the one-loop diver- 
gences in Weyl quantum gravity has been performed by 
Antoniadis, Mazur and Mottola [17] using methods de- 
veloped in [10]. The correctness of the /^-functions for 
the coefficients A and r\ calculated in [10] has been con- 
firmed. At the same time, the paper in Ref. [17] did not 
meet the suspicious J ^/^gR 2 counterterm. Indeed, this 
result coincides with our general expectations discussed 
above, but the situation with the two conflicting results 
does not look acceptable. In what follows we shall per- 
form a more general quantum calculation using dimen- 
sional regularization and starting from the action (2). In 
this way we will be able, in particular, to check the previ- 
ous calculations [10] and [17]. Even more important may 
be that we shall achieve better understanding of the role 
of the Gauss-Bonnet term in quantum gravity in n = 4 
and n = 4 — e dimensions. 

As far as one of our objectives is to perform a very 
complicated calculation in conformal quantum gravity, 
we have to provide maximal safety with respect to pos- 
sible calculational errors. For this end we shall use a 
new way of organizing calculations, which guarantees an 
efficient automatic verification of our result. Simultane- 
ously, we shall resolve another old-standing problem of 
quantum gravity. In the well-known paper [18], Capper 
and Kimber noticed that the Gauss-Bonnet term may, 
in principle, play a significant role in quantum gravity. 
Usually this term is disregarded because it is topologi- 
cal and does not affect the classical equations of motion. 
However, this conclusion is true only if the Bianchi iden- 



3 



tity is satisfied. This implies the diffcomorphism invari- 
ance of the theory. However, when the theory is quan- 
tized through the Faddeev-Popov procedure, the diffco- 
morphism invariance is broken and the vector space ex- 
tends beyond the physical degrees of freedom. In other 
words, after quantization not only the spin-2, but also the 
spin-1 and spin-0 components of the quantum metric be- 
come relevant, and the topological term may produce new 
vertices of interaction between these components. As a 
result, the quantum-gravitational loops may be, in prin- 
ciple, affected by the presence of the topological term. 
Of course, this output does not look probable, because if 
we include the topological term into the classical action, 
the gauge-fixing condition should modify and eventually 
compensate the new vertices. But this is a sort of be- 
lieves which are always good to verify. Such verification 
is one of the purposes of the present paper. 

We shall perform the one-loop calculation starting 
from the full action (2) , taking the topological term into 
account. As it was predicted in [18], the contributions 
of this term penetrate all vertices or, in other words, 
all elements of the background field method technique. 
However, despite many intermediate formulas strongly 
depend on the coefficient r\ , this dependence completely 
disappears in the final expression for the divergent part 
of the effective action in the n — > 4 limit. This cancella- 
tion means negative answer to the hypothesis raised by 
Capper and Kimbcr [18]. Moreover, it provides a very 
strong test for the correctness of the calculation. On the 
other hand, the 77-dependence is present in the ex- 
pression. Therefore, derivation of the relevant part of the 
effective action in an arbitrary dimension n^4 opens 
the way for constructing the complete 4 — e rcnormaliza- 
tion group equations in the conformal quantum gravity 
theory (2). As it will be shown below, the quantum ef- 
fect of the Gauss-Bonnet term leads to new fixed points, 
which have no analogs in the n = 4 case. 

The paper is organized as follows. In the next sec- 
tion we shall briefly describe the Lagrange quantization 
of theory (2). One can find a detailed description of this 
subject in Refs. [2, 10, 19, 20]. In section III the details 
of the bilinear expansion of higher-derivative gravity are 
presented. Some of the bulky expressions corresponding 
to this section are collected in Appendix A. Our expan- 
sions are more general than the ones which were known 
before [2, 11, 20], because they are performed for all 
higher derivative terms including ^J~gE and without 
taking into account the conformal gauge fixing condi- 
tion. This enables one, in principle, to derive divergences 
not only for the conformal case, but also for the general 
higher derivative quantum gravity [10-12]. In the present 
paper we perform only the calculation for the Weyl the- 
ory and expect to report the results for the general case 
later on. In section IV we derive the coincidence limit 
of the a 2 {x,x') coefficient of the Schwinger-DeWitt ex- 
pansion. The expression is obtained for the general n- 
dimensional space-time, in order to see the effect of the 
Gauss-Bonnet term more explicitly. After that, we de- 



rive the divergences of the Weyl gravity at n — > 4 and 
establish their independence on the parameter 77. In sec- 
tion V the renormalization group in the 4 — e dimensions 
is considered, and a number of new UV-stable and UV- 
unstable fixed points (due to the quantum effects of the 
topological term) are described. In the course of the cal- 
culations in sections IV, V we use the computer algebra 
program MAPLE (see, e.g. [21]). Finally, in the last 
section we draw our conclusions and discuss the possi- 
ble form of the non-conformal finite contributions to the 
one-loop effective action. 

II. QUANTIZATION AND GAUGE-FIXING 
DEPENDENCE 

The quantum gravity calculation in the background 
field method (see, e.g. [2] for the introduction) implies 
the special parametrization of the metric 

= 9\a> + h ^ ■ (6) 

In the r.h.s. of the last formula g^ v is the background 
metric and h^ u is the quantum field (integration variable 
in the path integral) . The I-loop contribution to the 
effective action of quantum gravity is defined as follows 
[10] 

= \ ^ DctW - l - In Det Y a * - 

- i In DetH gh , (7) 

where H is the bilinear (in quantum fields) form of the 
action (2) together with the gauge fixing term 

Sgf = M"- 4 / d n x^—g Xa Y a ? Xf3 . (8) 

The operator H g h is a bilinear form of the action of the 
Faddeev-Popov ghosts and \x is the dimensional constant 
(renormalization parameter in the dimensional rcgular- 
ization). The expression (7) includes also In Det V"' 3 , 
where V a/3 is the weight function. In the case of the 
higher derivative gravity theory this term gives relevant 
contribution to the effective action, because Y a/3 is a 
second order differential operator [10]. 

Introducing the gauge fixing term (8) one is fixing the 
diffcomorphism invariance. However, in the theory un- 
der consideration this is not sufficient, because there is 
another classical symmetry - local conformal invariance, 
which leads to a degeneracy even after the term (8) is 
introduced. Hence one has to choose the second gauge 
fixing condition. Following Fradkin and Tseytlin [10], we 
fix the conformal symmetry by imposing the constraint 
h = /ijJJ = 0. The interference between the two gauge fix- 
ing conditions may take place because the term (8) breaks 
the conformal symmetry in the background fields sector 
[16]. However, this breaking can not lead to the non- 
conformal countcrterms, because the last can be shown 
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insensitive to the choice of the gauge fixing condition. 
The general gauge fixing condition (here we restrict our 
attention to the linear background gauges) has the form 

V = \ {a^ a + tV^v,, - 6V v v„ + 

+ PiRfiv +P2-Rsv) , (9) 

where a, (3, 7, 5, p\ 1 P2 are arbitrary parameters. The 
action of the Faddeev-Popov ghosts has the form 



S gh = J d^x^f—g & (H gh )l C v , 



(10) 



where 



H gh = (H gh )l = -S;n- WW, 2/JV M V . (11) 

The parameter (3 is fixed in the conformal case due to the 
conformal symmetry condition = 0, hence (3 = — 1/n. 
Other parameters may take different values and their 
choice may influence, in principle, the one-loop diver- 
gences. The general analysis [16] shows that the C 2 , E 
and R 2 counterterms can not depend on these parame- 
ters while the D-R-type counterterm may have such de- 
pendence. In what follows we shall use these data exten- 
sively, namely we will not pay attention to the irrelevant 
J yj—g OR counterterm and, on the other hand, we shall 
choose the gauge fixing parameters a, 7, 5, pi, P2 such 
that the calculation of other counterterms becomes sim- 
pler. Let us remark that the dependence on the param- 
eters pi, P2 has been explored and found irrelevant in 
[12] for the non-conformal version of the higher deriva- 
tive quantum gravity. 



III. BILINEAR EXPANSION 

The action (2) includes only higher derivative 
conformal invariant and surface terms. There are no 
Einstein-Hubert, cosmological and J \J~ r g R 2 terms in 
the action, because none of them possesses local confor- 
mal symmetry. But for the sake of completeness, the 
bilinear expansions for all these terms will be given too. 
The parametrization of the quantum metric h, v will be 
chosen according to (6). Let us remark that the rele- 
vant divergences in the theory (2) are independent of the 
choice of parametrization for the quantum metric [16]. 
When making the expansions of the elements of the grav- 
itational action, we keep in mind that the relevant terms 
are of the second order in h, v . Hence we shall pay the 
main attention to this order of the expansion. In what 
follows we indicate all quantities constructed from the to- 
tal metric g' by prime (e.g. g'^ v 



etc), and reserve simpler notations (e.g. g 1 *", \f—g, r^„, 
R a ,p v etc) for the quantities constructed from the back- 
ground metric g, u . 

For g'^ v and ^j— g' the expansions can be presented 

as 



tT - + h^h Xu - h^h x h TV + ... 
% + .9(1) + .9(2) + .9(3) + - 



(12) 



and 



-.9 (l + \h+ l -h 2 -\h^ + ...). (13) 

For the coefficients of the affine connection, using (12), 
we arrive at the expansion 



C - + ST» U , where 51%, = £ ST™ (14) 

n=l 

Here the tensors ST^ ™ are given by the expressions 

*T (n) % = \ 9(n-i) (v„ V + V, V - V/3 V) • (15) 

For the curvature tensor one can establish the following 
general expression: 

00 

- sv^ sr^ = + E R{n) V ■ ( 16 ) 

n=l 

In the first and second orders in the quantum metric, 
, we obtain the following expressions for the Riemann 
tensor: 

- Vn^hpv + R\, v hp - R X -,, U ft-") 1 

R {2) %„ = \h aX (v t y x h ul3 + V u Vph ll x + 

+ VvVphxp) + \ [v^« A (V A /^ - 

- Vph vX - V„^ A ) + Vphi (Vxh% - 

- V a h Xfl ) + V v h x (Vxh^- W a h x ,) + 

+ V x h^(VxK- V a ft A ,)] -{n~v). 

(17) 

For the Ricci tensor, similar expansions have the form 
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+ i v Q v (v Q h£ - v"*c) + i v u h afi + 1 (v^ + v, v - v) • ( 18 ) 



For the scalar curvature we meet the following expan- 
sions: 

R {1) = V^VuhT - oh - h^R^, 



R( 2) = h a ( , (V a Vf,h + nh a p-V a Vxh%-V x V a h%) + 

1 
4 



+ V a h aX (V x h-Vph l3 x )-- A V x hV x h + 



h aX h x R a ^ + -V x h aP V x h 

\v a h X[) v p h aX . 



a/3 



(19) 



With these expansions in hands, we can derive the part 
of the action quadratic in the quantum fields. It proves 
useful to consider an alternative version of the action (2) 

S w (n) = V"- 4 > J d n x^7j {xRl vaP + yRl„ + 

+ zR 2 + tDR} , (20) 

where the new parameters x, y and z are related to 77 
and A as follows: 



2A 



(n — 2) A 



4 77, 



z = 1] + 



1 



A(n-l)(n-2) 
After some algebra we arrive at the formula 



(21) 



S (2) 



-n 



(n-4) 



(V~9 R%) + z (y=g R 2 ) } , (22) 



- D 2 \ (2) 

9 R^a/3 J + 

(2) 



y 



where the complicated expressions for the bilinear forms 
are collected in the Appendix A. 

Starting from the expression (22), and using (A4),(A5) 
and (A6), one can easily find the bilinear form of the 
action (20). The operator Ti. depends on the gauge fixing 
term (8). The gauge fixing parameters a, (3, 7, 5, pi, P2 
in (9) will be chosen in such a way that the operator takes 
the most simple, minimal form 



where K is a non-degenerate c-number operator. The 
two of the possible non-minimal fourth derivative terms 
g^vVaOV/3 and (/q^V^DV^ vanish due to the conformal 
gauge fixing condition hit = 0. The simplest choice of the 
parameters providing the cancellation of the remaining 
non-minimal fourth-derivative structures V^V„V Q V/3 
and g^V^DVa is the following: 



a 



y + 4x 
1 , Pi 



1 

P2 



2x-2z 
y + Ax 
0, 



(24) 



with /3 = — 1/n defined by the conformal gauge fixing (as 
we already noticed). Let us remark that this "minimal" 
choice of the gauge fixing is sensitive to the introduction 
of the Gauss-Bonnet term, as we expected. On the other 
hand, if we fix the value of 77 such that the sum of the 
Wcyl term and the topological term gives 



C 2 - E = 2W = 2(i£„-ifl 2 ), 



(25) 



the gauge fixing condition (24) coincides with the one of 
[10, 17]. 

After some algebraic calculations and using the com- 
mutators (A3), we find 

I s + s af] (2) = h ^ H hafi > where 



H 



KU 2 + D px V p Vx + JV"V„ - (V M Z") + W 



(26) 



and K, D pX , N^V^, V^Z^ and W are matrices in the 
/V^-spacc. In order to derive an explicit form of A M 
and one has to extend the derivation of bilinear ex- 
pressions of the Appendix A. However, the derivation of 
these quantities does not have much sense, because the 
terms iV M V M and (V^Z^) may be safely disregarded. 
The reason is that both expressions and Z^ are 
covariant derivatives of curvatures. Therefore they may 
contribute only to the irrelevant gauge-fixing dependent 
/ ^/^gOR-type counterterm, which we are not calculat- 
ing here. Below we shall simply set both terms to zero. 
Let us introduce the useful notation 



H = Kn 2 +o(v 2 ). 



(23) 
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for the projection operator into the traceless sector of the matrix. Without this condition the identity matrix is 

W space. Here, as usual, S pua p. Taking the conformal gauge into account, we 

^ find 

fi^u.ap = 2 (sw g»p + g^p gua) ■ 

{K)^, a f} = (| +z) V,a/3, (27) 

Since we assume the conformal gauge fixing condition 
h = 0, the tensor <5 M „ iCe( 3 plays the role of the identity 



(D px )^ afj = -2xg uf} R a f\ + AxS p R x a ^ + (3s + y)g px Rp tav p + 2x5„p pX R fia - {Ax + 2y)8 p a R x „g u/3 - 
- 2xg v(i R iia g px + ^^5^ a0 R pX - zg^5^ a px R + ^z5^ a[j g px R - 2z8 a[j px R^ ; 



and 



iW)^v,ap = -^-gvf3R t j.pXaRa PX ' 7 H ^ R f ' ap, Ru/3pX H 7, R ap P Rxvfip H 7, R P ' vRpaXfi + 



3x 
~2 

y — 5x 



- — ^ RppR P au(3 + RpaRvfi + -^gu[)RppR' 'a + -^-gv[jRRap — -jRRvfiap + 



+ zR^R af3 - - (xR 2 pXar + yR 2 pX + zR 2 ) (8^ ta0 ) . 

I 



(28) 



(29) 



In the above formulas we used special condensed nota- 
tions which enable one to present the expressions in a 
relatively compact way. The idea of these condensed no- 
tations is that all the algebraic symmetries are implicit, 
including the symmetrizations in the couples of indices 
(a/3) <-» {fiv), (a <-» (3) and (//<-» v) , and also in the 
couple (p <-» A) in the operator D pX . In order to obtain 
the complete formula explicitly, one has to restore all the 
symmetries. For example, 

RfipR P aisfi * ^ (RppR^ 1 av(3 ~t~ RapR P p*fiv) 

restores the (a/3) <-> (/if) symmetry The same pro- 
cedure has to be applied also for the other symmetries 
(p-*-»A), (a <-> /3) and {p^v). 

In order to use the Schwinger-DeWitt method for the 
fourth derivative operator [10], we need to reduce it to 



the minimal form (23). For this end one has to multiply 
the operator (26) by the inverse matrix K^ 1 , given by 



1 \pv . a/3 



y + Ax 



fip,v , a/3 



Let us notice that the matrix K 1 is a c-number op- 
erator and hence this multiplication does not affect the 
divergences. By straightforward algebra, one can find the 
minimal operator 

H = R- 1 H = ID 2 + V px \7 p \7 x + U , (30) 

where the new expressions 

ypx = K-^t) pX and U = K~ X W (31) 
already do not possess the symmetry in (a/3) <-> {[iv). 
The expressions for these two matrices are the following: 



\y)nv,aP — y _|_ 4^ 1 ~~^~9vP u p,p\<T U a H 2 a/i u \v0p-\ ^ p » u paXp ^ ^ H PP U ctv0 + 



+ — Rp a R v p + -^gvpR x ^R a \ - 4 (xR 2 pXaT + yR 2 prJ + zR 2 ) (Spv , a/3 ) + 

+ -^gvflRRap H ^ -^ P «M RvfipX + zRpvRafi ~ -^RRvPctp, \ , (32) 



4 where 



T> A = — V>iki, (33) 

y + 4x ,=i ^=5^^^; k 2 = V, Q/3 ^ A i?; (34) 
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= &vti pX Rg^a] k 6 = 5^ , a/3 R pX ; 
k 7 = i(^i? A ) a/3Al + ^i? A )^ a ); 
ks = S^^-R^a) 5 k 9 = g v/ j i?( Q pA Al) ; 



kio — -^(S a p p R^v + 5^u P Rap) 



and 



61 = -2x; 6 2 = z/2; fe 3 = 3x + y; 

64 = 2a; 65 = — z be = x + y/2 ; 

67 = —4a; ; 6g = — 4x — 2y ; fog = —2a; ; 

610 = -2z. (35) 

The above form of V rpA is helpful in organizing the cum- 
bersome calculations of divergences which will be de- 
scribed in the next section. 



IV. DERIVATION OF DIVERGENCES 

The algorithm for the 1-loop divergences correspond- 
ing to the minimal fourth order operator can be written 
as [10] (here we use the Euclidean signature of the metric, 



in order to be consistent with [10]) 



- In Det H 
2 



where 

lim a 2 (x' , x) 



1-4 



x 



div (47r) 2 (n-4) 

x / d n Xy/g tr lim a 2 (x',x), (36) 

J x'^x 



— R 2 

90 c""! 3 



+ -/v,,/v'"" U • j^A'V ' (37) 



l R ^ PX + ^V p p V\ + lv pX y P X. 



Here TZ^ is the commutator of the covariant derivatives 
acting in the tensor h a/3 space, 



% v = [V M , V,] . 



(38) 



The full collection of the traces of the expressions (37) is 
presented in the Appendix B. 

It is straightforward to find the contributions of the 
weight and ghost operators (the algorithm for the non- 
minimal vector operator can be found in [10, 12]) 



-- In detF| d i- 



—i In detHgh |div = - 
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(n-4)(47r) 2 



(39) 



(n-4)(4 7 



1 ,2 1 , 17 

- £ f H 

6 s 3 ^ 30 



)*}■ 



(40) 



where the parameter £ is given by 

n- 2 



2(n- 1) 



(41) 



Collecting all the results from (37), (39) and (40) ac- 
cording to (7) and using the formulas from Appendix B, 
we arrive at the functional trace of the overall coincidence 
limit of the a 2 (x', a;)-cocfficient 



4* 



lim sTra 2 (x',x) = lim Tr a 2 (x' , x){H)- 

x'^x x' —*X . 

- Tra 2 (x\x)(Y)-2T±a 2 (x',x){H gh )\ • (42) 



The last expression can be regarded as a functional super- 
trace of the coincidence limit of the a 2 (x', a;)-coefficicnt 
of the differential operator acting in the direct product 
of the tensor h^, vector (third ghost) and vector ghost 



spaces. The sign difference between the different terms in 
(42) is due to the different Grassmann parity of the fields, 
and the operator Tr includes integration, as usual. 

Let us present the result in terms of the parameters r\ 
and A : 

A\ = -fj, n - 4 j ' d n xj—g {p 1 (n)E+ 

+ f3 2 (n)C 2 +(3 3 (n)R 2 } , (43) 

where the coefficients (/3- functions) /3i(n), (3 2 (n) and 
/3s(n) are given by the expressions 

p t (n) = 6^ +6<pr, + 6? ) r ?, < = (1,2,3). (44) 

The coefficients 5^ are the following functions: 
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6<f 
4° 



si 2 



AO 



A 2 



15 n 6 + 86 n 5 + 201 n 4 - 4842 n 3 + 8104 n 2 + 6624 n - 9648 



2880 (n- 1) (n- 3) 2 



(n - 4) (n - 2) (r 



39n 3 -40n 2 -196n+192)A 



48n(n- 1) (n - 3) 2 

(n-4) 2 (n 3 + 9n 2 + 14n+12) (n - 2) 2 A 2 
48 (n - 3) 2 n 2 

i-2) (5n 6 + 299 n 5 - 1162 n 4 - 2570 n 3 + 15056 n 2 - 18528 n + 6720) 



960 n(n- l)(n-3) z 
n - 4) (n 4 - 3 n 3 + 50 n - 36) (n - 2) 2 A 



48n(n- 1) (n - 3f 
4) (n 3 + 10 n 2 + 10 n + 24) (n - 2) 3 A 2 



48 (n-3) n 2 
n - 4) (5 n 5 + 22 n 4 + 179 n 3 - 930 n 2 - 112 n + 816) 



960 (n - l) 2 (n - 3) 
4) (n 4 - 4 n 3 - n 2 + 10 n - 12) (n - 2) 2 A 



24n(n-iy (n-3) 
n - 4) 2 (n + 1) (n 2 + 2 n + 12) (n - 2) 3 A 2 



96 (n - 1) (n - 3) n 2 



(45) 



r 



The above coefficients, despite their chaotic appear- 
ance, provide a lot of important information. First of all, 
they show that for the Gauss-Bonnet topological 

term contributes to the effective action in a non-trivial 
way, and in particular produces the / y/^gR 2 -type term. 

On the other hand, it is remarkable that all 5^ and 
<5- 2 ^ coefficients are proportional to (n — 4). Hence, for 
n — 4 one can see that the n-dependence completely 
disappears, and the final result, Eq. (43) becomes very 
simple. Let us write down the expression for the one-loop 
divergences 



div 



M"~ 4 f ,n , — f 137 ^ 199 _„ 



(47r) 2 (n-4) 



y \ 20 30 



(46) 



where we used Eq. (25) and the pseudo-Euclidean signa- 
ture. The expression (46) coincides with the one derived 
by Antoniadis, Mazur and Mottola in [17]. Both coeffi- 
cients in (46) also coincide with the ones derived by Frad- 
kin and Tseytlin in [10]. However, we do not meet the 
polemical / ^/^gR 2 -type divergence [10] and hence there 
is no need to apply the conformal regularization [13, 14] 
discussed in the Introduction. As far as our calculation 
is seriously tested by the cancellation of the numerous n- 



dependent terms in the n — > 4 limit, we strongly believe 
in its correctness. Thus, the / y/^gR 2 -type one-loop di- 
vergence does not show up in the one-loop effective action 
of Wcyl quantum gravity. 

The expression (46) does not contain the divergences 
in the cosmological constant and of the Einstein- Hilbcrt 
type. This fact is due to our choice of the regularization 
procedure. As far as we start from the conformally invari- 
ant action, the theory does not possess dimensional pa- 
rameters and therefore the divergences of these two sorts 
may be only quartic and quadratic ones. Of course, in 
the dimensional regularization which we are using here, 
the quartic and quadratic divergences do not show up. 
However, one can easily see the cosmological and linear 
in R divergences in other regularization schemes, for ex- 
ample in the covariant cut-off [10] or in the covariant 
Pauli-Villars [8] regularizations. Of course, the logarith- 
mic divergences which we have calculated (46) do not 
depend on the choice of the regularization scheme. 

Including the matter fields one meets additional con- 
tributions to the divergent coefficients in (46). As it 
was already noticed in the Introduction, the conventional 
scalars, fcrmions and vectors give contributions of the 
same sign to both f3\ and /?2 , while the contributions of 
higher derivative scalar and fermion have opposite sign 
[4-6]. The sign of the coefficients in (46) coincides with 
the one of the conventional fields. Hence, since the J R 2 - 
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type divergence is absent, one can use the method of [5] 
and adjust the number of higher derivative scalars and 
fermions such that the one-loop divergences completely 
cancel. In a more complicated situation, when the mat- 
ter coupled to quantum gravity possesses self-interaction, 
the quantum gravitational effects modify the divergences 
and the renormalization group trajectories, also, in the 
matter sector of the theory. This issue has been studied 
in details [15, 23, 24] for the case of the higher derivative 
gravity [33]. In both conformal [15] and general [23, 24] 
cases the effect of quantum gravity is rather smooth 
and always favors the asymptotic freedom in the mat- 
ter fields sector. Due to the absence of the / R 2 -type 
divergence the investigation for the conformal case can 
be indeed performed without the special conformal regu- 
larization (which we discussed in the Introduction), while 
the quantum-gravitational corrections to the /3-functions 
in the matter sector are exactly the ones derived in [15] 
(see also [2]). The reason is that these (3- functions do 
not depend on the scalar curvature and hence arc not 
affected by the conformal regularization. 



V. RENORMALIZATION GROUP EQUATIONS 

The renormalization group (RG) equations for the the- 
ory (2) may be considered in two different ways [34]. 
The first possibility is to take usual n = 4 /3-functions, 
in this case we meet exactly the same RG equations as 
in [10]. It proves useful to introduce a new parameter 
p = —1/rj. Let us remark that the choice of A as a 
coupling constant in the action (2) is fixed, because (i) 
A is a parameter of the loop expansion in this theory; 
(ii) One can not change the sign of A without chang- 
ing the positivity of the graviton energy. At the same 
time there are no similar constraints for the coefficient 
of the Gauss-Bonnet term and therefore the choice can 
be made according to convenience. The usual n = 4 
renormalization group equations for A and p have the 
form 



d\ 


d\ 




~dt 


^ dp 




dp 


dp 
^ dp 




~dl 




where 







&(4) - -a 2 X 2 , A(0) = A ; 



= 0,(4) = -b 2 p 2 , p(0)= Po , (47) 



199 



15 (4tt) 2 



b 2 = 



261 



60 (4tt) 2 



(48) 



The above equations indicate the UV asymptotic freedom 
in both parameters. In other words, there is a single fixed 
point A = p = and it is stable in the high energy limit 
t — > oo. 

Let us now consider a more complicated version of the 
renormalization group equations, taking the dimension 
n = 4 — e for — 1 < e < 1. Mathematically this means 
that we do not take the limit n — > 4 in the equations 



(47) . The renormalization group equations which emerge 
as a result of this procedure will be different from (47) 
and one can expect to see qualitative effects of the Gauss- 
Bonnet term in this framework. 

Similar approach to the renormalization group proved 
fruitful in the two-dimensional quantum gravity [25] , due 
to its relation to the concept of asymptotic safety [26] and 
to the discussion of the universality classes of quantum 
gravity theories [27]. The main idea of 2 — e quantum 
gravity is the following. In the precisely n = 2 dimen- 
sions, quantum gravity is a topological theory similar to 
the one which we meet in n — 4 starting from the pure 
Gauss-Bonnet term. But, if we generalize the theory for 
n = 2 — e , there is a dynamics (different from the Gauss- 
Bonnet theory, where the propagator does not appear 
even for n ^ 4) and at the quantum level one meets a 
non-trivial UV fixed point of the renormalization group 
[25-27]. Keeping this example in mind, it is natural 
to expect that the effect of the Gauss-Bonnet term on 
the renormalization group equations in n = 4 — e may 
be non-trivial and in particular may produce new fixed 
points. 

Consider the renormalization group equations for A 
and p in n = 4 — e dimension. The naive form of the /?- 
functions would be based on the "standard" expressions 
(47) 



-eA + /3 A (4) : 



A. 



-ep + (3 p (4), (49) 



indicating the one extra nonzero fixed point for each of 
the effective charges A(t) and p(t). Indeed, the fixed 
point A = p = remains stable in UV for e > 0. How- 
ever, this naive consideration is incorrect because the 
Gauss-Bonnet term gets dynamical in n ^ 4, affecting 
the renormalization group equations in a non-trivial way. 
Using the expressions (45), we arrive at the following cor- 
rect form of the renormalization group equations, quite 
different from (49): 

= ^ A "(4rfl^% +ff V)- (50) 



dX 
~dl 



The coefficients f 1,2,3, .91,2,3 may be expressed via the 
coefficients from (45) as 



c(0) 



f2 = 5{ 1) /\, f 3 = 6?>/\\ 



(2) 



91 



= 4 0) , ff 2 = 4 1) A, <?3 = 4 2) /A 2 . (51) 



One can remark that f 1,2,3 and 31,2,3 depend only on 
the parameter e and not on the couplings. In the limit 
e = we come back to the equations (47). The renor- 
malization group equations (50) are non-linear and do 
not admit a simple analytic solution. For this reason we 
shall start from the search of the fixed points, that are 
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the values of A and p for which both /3-functions van- 
ish. Consequently, we explore the stability of these fixed 
points and establish the renormalization group flows for 
some particular values of e. 

In order to find fixed points, we consider the particular 
values of the parameter, e = 0.9 , e = 0.1 , e = 0.01, 
e = —0.01 , e = —0.1 and e = —1. The numerical values 
of the coefficients for these cases are presented in the 
Table 1. The point e = 0.9 is numerically close to e = 1 
(n = 3) where the expressions for the /3-functions become 
singular. 



e 


h 


h 


h 


9i 


92 


53 


0.9 


-16.77 


-28.73 


-36.48 


2.359 


-42.51 


-46.98 


0.1 


-4.301 


0.08 


-0.016 


6.385 


-0.174 


-0.318 


0.01 


-4.344 


0.008 


-0.0001 


6.608 


-0.016 


-0.03 


-0.01 


-4.356 


-0.008 


-0.0001 


6.659 


0.016 


0.03 


-0.1 


-4.416 


-0.086 


-0.013 


6.902 


0.146 


0.286 


-1 


-5.416 


-0.947 


-0.81 


9.98 


1.087 


2.526 



Table 1. Numerical values of the coeffi- 
cients for the particular values of e. 

The numerical analysis shows that for each of the 
choices e = 0.9, e = 0.1 and e = 0.01 there are four 
fixed points which are new compared to the e — case; 
while for the values e = —0.01 , —0.1 , —1, there are two 
new fixed points. The values of the parameters corre- 
sponding to these fixed points are shown in the Table 
2. 



Fixed Point for e = 0.9 


1 


2 


3 


4 


Ai 





-6.817 


-5.945 


1.807 


Pi 


-8.475 


-10.75 


-12.52 


-3.05 


Fixed Point for e = 0.1 


1 


2 


3 


4 


A l 





-14.421 


-1.232 


16.236 


Pi 


-3.671 


-3.159 


-3.647 


-3.706 


Fixed Point for e = 0.01 


1 


2 


3 


4 


Ai 





-5.228 


-0.119 


5.457 


Pi 


-0.364 


-0.351 


-0.363 


-0.371 


Fixed Point for e = -0.01 


1 


2 


3 


4 


Ai 









0.1186 


Pi 


0.3625 






0.3628 


Fixed Point for e = —0.1 


1 


2 


3 


4 


Ai 









1.147 


Pi 


3.576 






3.597 


Fixed Point for e = — 1 


1 


2 


3 


4 


Ai 









8.001 


Pi 


29.157 






30.239 



Table 2. Numerical values of the param- 
eters corresponding to the new fixed points. 
No one of them has analog in the e = case. 

The stability properties of the fixed points can be eas- 
ily investigated in the linear approximation. The result 



is that, in the cases e = 0.9, e = 0.1 and e = 0.01, the 
fixed points (Ai,pi) and (\2,P2) are saddle points while 
the fixed points (A3,p3) and (X4, p^) are absolutely un- 
stable in the UV limit t — > 00. It is worth noticing the 
contrast with the e = renormalization group equations 
(47) with a single UV stable fixed point A = p = 0. In 
the n > 4 cases e = —0.01 , e = —0.1 and e = — 1 there 
are only two extra fixed points, one of them UV-stable 
and IR-unstable and another one a saddle point (unsta- 
ble in both UV and IR regimes) . As shown in the Figures 
1 and 2, there are no additional (compared to the stan- 
dard e = case) UV-stable fixed points for positive e, 
and at the same time, for negative e there is always one 
additional fixed point with stability in the UV domain. 

An interesting observation concerning the renormaliza- 
tion group trajectories is that no one of them crosses the 
line A = 0. In other words, the renormalization group 
flow in this theory is divided into two separate parts: one 
with A > corresponds to the positively defined energy 
of the gravitons (massless spin-2 mode) and another one 
to the unphysical graviton sector with A < 0. There are 
examples of the qualitatively new UV-stable fixed points 
with A > (see Figures 3,4). At the same time there 
are no such examples for the case A < 0. One can sup- 
pose that this property of the fixed points is related to 
the limit e — > 0, where all new UV-stable fixed points 
presumably should tend to A = 0. 

It is obvious that no one of the fixed points which we 
have found so far coincides with the standard one A = 
p = of the n = 4 renormalization group. The natural 
question is whether it is true that the effect of the Gauss- 
Bonnet term is to kill the asymptotic freedom in n = 4— e 
dimensions. The answer to this question is definitely not. 
The source of our failure to see the standard fixed point 
is that we have used only the algebraic equations (3\ — 
(3 p = and due to the non-polynomial form of the (3\- 
function (50) one can not see the fixed point with p = 
in this way. So, in order to complete our study we have to 
consider, especially, the possibility of simultaneous A — > 
and p — ► . Using elementary transformations, one 
can check that the regimes A</> and p«A lead to 
contradictions. Therefore, we consider, additionally, the 
possibility of the special solution p = kX, where k is a 
constant. Under this assumption, the equations (50) are 
consistent if 

h k 3 + (2<?i - h) k 2 + (f 3 - 2g 2 ) k + 2g 3 = , (52) 

with an additional restriction 
2 



b = 



(4irk) : 



(gik 2 -g 2 k + g 3 ) > 0, 



(53) 



dictated by the asymptotic freedom, in the UV regime. 
The origin of this condition is the following. After the 
relation p = k-X is imposed, the equation for A becomes 



dX 
~dt 



= -eX - bX z 



(54) 
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FIG. 1: Diagrams for e = 0.9. The j/-axis represents the coupling A and the rr-axis, the coupling p, as well as in all subsequent 
plottings. The left diagram shows the fixed points 1 and 4, and the other shows the points 2 and 3. The labels of the fixed 
points correspond to the numeration in Table 2. The arrows indicate the direction of the renormalization group trajectory at 
the given point (we have also drawn some trajectories for illustrative purpose). One can distinguish stable, completely unstable 
and saddle fixed points at these and further diagrams. 




FIG. 2: The fixed points 1 and 3 are shown for the cases e = 0.1 (left) and e — 0.01 (right). Clearly, the point 1 is a saddle 
point (unstable) and the point 3 is UV-unstable. The points 2 and 4 are similar (saddle and UV-unstable, respectively) and 
are not plotted. 



The general solution of this equation has the form 

X{t) = b{e«- 1) ' A ^)>°- ( 55 ) 

It is easy to see that the asymptotic freedom in the UV 
limit t — > +oo holds for e > and (53) satisfied. For 
e < and condition (53) satisfied the situation is more 
complicated, because the UV stable fixed point is non- 
zero X(t — > oo ) = — e/b > 0. However, this fixed point 
tends to zero when e — > 0, and we can consider that the 
theory is asymptotically free in this sense. At the same 
time, independent on the sign of e, the theory with b < 
does not manifest the asymptotically free behavior in the 
UV limit. 



In fact, there is no guarantee that the condition (53) 
is satisfied for every choice of e and all real roots of 
(52). The numerical computations show that for e = 
—0.01 , —0.1 and —1, the equation (52) has only one 
real root and that this root satisfies the condition (53). 
On the other side, for e = 0.9 , 0.1 and 0.01, the equation 

(52) has three distinct real roots, one of them violating 
condition (53) for each case. Thus, there are solutions of 
the equations p = kX and (52) which do not satisfy Eq. 

(53) . Let us remark that, in all cases we have examined, 
there are also solutions with the UV stable fixed point 
(0, 0). However, the asymptotic freedom depends on the 
choice of the initial condition on the A— p plane. In some 
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FIG. 3: Fixed points 1 and 4 for the cases e = —0.01 (left) and e = —0.1 (right). The point 1 is a saddle point (unstable) and 
the point 4 is UV-stable, contrary to the analogous points with positive e. 




FIG. 4: The case e = — 1, with a saddle point (1) and a 
UV-stable one (4). 

cases, when Eq. (53) does not hold on the special solu- 
tion of the equations (50), the (0, 0)-point is not stable 
in the UV. 

Looking at the Figures 1-4, one can observe the renor- 
malization group trajectories (for the e > case) linking 
the IR-stable point (A3, pz) to the UV-stable point (0, 0), 
or alternatively the IR-stable point (A4, pi) to (0, 0). The 
situation is similar for e < 0, but here the renormaliza- 
tion group flow is inverted, linking the IR-stable point 
(0, 0) to the UV-stable point (A4, p±). 



VI. CONCLUSIONS AND DISCUSSIONS 

We have calculated the one-loop effective action for the 
Weyl gravity with the Gauss-Bonnet term. In the n — > 4 
limit the quantum effects of the Gauss-Bonnet term can- 
cel. This cancellation may be seen as a negative answer 



to the problem raised in [18]. This result is valid, at 
least, in the framework of the conformal quantum grav- 
ity. Another remarkable fact is that, in agreement with 
[17], there is no infinite J yJ—gR 2 counterterm. Other 
sectors of the divergent part of the effective action are in a 
perfect agreement with both earlier calculations [10, 17]. 

Despite the one-loop divergences are conformal invari- 
ant, this symmetry may be broken at the one-loop level in 
the finite part of the effective action. The divergences of 
the J \J—gC 2 and / y/—gE-type produce the anomalous 
violation of the Noether identity (4), and as a result the 
finite part of the one-loop effective action contains usual 
non-local [35] anomaly-induced terms [4]. There may be 
also local J s/^g~R 2 -type contribution which deserves a 
special discussion. It is easy to see that there are two dif- 
ferent possible sources of this term in the Weyl quantum 
gravity: 

i) If the calculation is performed in a dimensional reg- 
ularization, the 6^ and 5^ terms in (45) are propor- 
tional to (n — 4) and therefore they produce the finite 
f ^J—gR 2 term directly from A\. It is remarkable that 
this contribution depends on the coefficient -q of the 
Gauss-Bonnet term. According to [8], this contribution is 
a subject of strong ambiguity typical for the dimensional 
regularization. In general, the dimensional regularization 
is unable to predict any definite value for the coefficient 
of the finite / yj— gR -term. 

ii) The infinite J y/—gDR-type counterterm, which we 
did not calculate here, may produce contribution to the 
conformal anomaly and eventually to the finite J ^J—gR 2 
term. However, this contribution is plagued by a double 
ambiguities. First of all, the J y/—gOR-type countert- 
erm itself is gauge-fixing dependent [16]. As it was al- 
ready explained above, this is the reason why we did not 
calculate this counterterm. The second source of ambi- 
guity is the derivation of anomaly and of the anomaly- 
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induced effective action. In relation to the J yJ—gR 2 
term these procedures may be ambiguous. The detailed 
discussions of this issue have been given recently in [8], 
where the ambiguity has been confirmed not only for 
the traditional version of the dimensional regularization 
(where it is completely out of control) but also in a more 
physical covariant Pauli-Villars regularization with non- 
minimal scalar massive regulators. It is worth noticing 
that the status of this last ambiguity in the Weyl quan- 
tum gravity is quite different from the one in the semiclas- 
sical approach. In the last case the ambiguity is always 
reduced to the freedom of adding the / \J—gR 2 term to 
the classical action of vacuum, while in the former case 
this operation would increase the number of physical de- 
grees of freedom (see, e.g. [2] and references therein) 
and hence can not be seen as the legal operation for the 
theory (2). 

In any case the local conformal invariance in Weyl 
gravity is violated at the one-loop level by quantum 
corrections. Hence, despite the general higher deriva- 
tive quantum gravity is indeed renormalizable [19, 29], 
the particular conformal version is multiplicatively non- 
renormalizablc at higher loops. Our results show, how- 
ever, that the conformal quantum gravity can be re- 
garded as a good approximation. The corresponding pro- 
cedure means that one can start from the theory with a 
very small coefficient of the / \J—gR 2 term. Due to 
the one-loop renormalizability of the conformal theory 
this coefficient will remain very small at the quantum 
level. If we consider the conformal quantum gravity in 
this framework, the problem of ambiguity of the anoma- 



lous J yf—gR 2 term is irrelevant and we can regard this 
theory as a useful particular example of the higher deriva- 
tive quantum gravity models. 

One of the outputs of our investigation are new fixed 
points of the renormalization group flows which appear 
due to quantum effects of the topological Gauss-Bonnet 
term in 4— e dimension [36] . One can expect even greater 
number of non-trivial fixed points for a general higher 
derivative quantum gravity, with the Einstein-Hilbcrt, 
cosmological and J i? 2 -terms included. 
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APPENDIX A: BILINEAR EXPANSIONS 
QUADRATIC IN CURVATURE TERMS 

In this Appendix we collect the cumbersome expres- 
sions with the bilinear expansions of the relevant terms 
of the second order in curvature. Furthermore we present 
the transformation of these terms to the form which is 
useful for the derivation of the effective action. The ini- 
tial set of the bilinear expansions has the following form: 



(2) 

(V=g R 2 pva p) = yf^hT { <5„„ , Q(3 n 2 - g pa V^aV„ + V a VpV p V„ - g lux V p V p V v V p + 4:R p av pV p V p + 
+ <V , af3 R P xV p V x + R palJ pU - 2R p ^ av VpV p - 2g^R paX pV x V p - 4g uf3 R apXf ,V x V p + 

+ -^9vuRn P \TRf3 PXT + 9 P vR p \TaR pXT — J {$nv,aP — 2 9n»9af3jR 2 p \ T 8 — 2 -^WvA-Ri//^ j t 

(y=9 R%) ^ = V=9 h pv { iv Q V p V fi V„ + \s P v, Q/3 n 2 + ^aVpV^V/jV - ^a/sVpV^V - 

- ^^VaV^V^V + iff^5a/3V p nV p - ^g av VpV p a + ^g pv V a VpU - ^g„ a n\7 p \7 p - 

- 2g va R p V p V p + -6 pv . a pR px V p V x + g a pR p „V p V p - R p pV a V„ + 

+ g va R pS 3^ + 2g va R p M V ' P V f3 + g pu R p pV i?" A V p V A - 2g p pR p „V a V p + 

+ g (g^ga/3 — 2<5 M „ , a p)R 2 p \ + RfiaRvf) + 2g va R^ p R f ' p — g a pR PP R f V j h af3 ; 



(Al) 



and 



(V^i? 2 ) = ^gh pv {VnVvVaVp-gapV^VvD-g^DVaVp+g^gapn 2 ~g va RVpV p - 
- 2R^X7 Q V p + g pu RS7 Q V p + 2g a pR^U + ^(S pViCt p - g tlv g a p)RO + 
+ lRgv(iR P a — g P vRRot0 + -^(g^g a p — 2<5 M „ , a p)R 2 + R pu R a p^ h afj . 



(A2) 
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It proves necessary to establish some commutation re- 
lations between covariant derivatives. In the expressions 
below we have omitted those terms which may contribute 
only to the total derivatives in the effective action. Also, 
for the sake of brevity we broke the symmetries in the 
pairs of indices (a/3) and (fJ,v). These symmetries has 
to be restored for practical calculations. We can write 

5w3V A V M V a V A h al3 = {g ul3 R pt y»V a + R^x^V x V a + 

+ ^aVX) h a ? , 

g v pUV a V p h aP = {gvpRncP + R^ afl a + 

+ g vl3 V ll uV a +g vl 3R pi y , V a + 
+ 2R ufjx ,V x V a ) h a ? , 



5^V P V A VV a h a P = (R pX V?V x + a 2 )h, 



+ 2R pjCt V v Vp — R p ' vap y ' pV '/3 



VcV^V^/i^ = [V iU V i ,V a V / 3+4 J R iUa V I/ V^ 
- i? p raAI (V (9 V (3 + V /3 V p ) + 
+ 2R" a ^V p V v ] h^, 



g,pV x V a V p V x h a P = {g„pR PP V p V a +R^x P V x V a + 
+ g v pV p O\7 a +g„ p R x a \7 x \7 p + 
+ Rv l 3axV x X7 p +g„ l3 R l , a n + 

+ Rvfian^ + -^gvf}RfiT P \R a Tpx + 

+ \R a / x Rvf ip x)h ati , 



g vP V a UV p h al3 = (2g„ fj R p aXp V P V x + 

+ g vfj (2R pa V p V p + R a „n + 
+ V^UV a + 2R l3uXa V p V x )h a \ 



V x V„V x V„h = (R PP X7P\7 lJ -R pljXp \7 x \7P + 

nVpVvh = (2R p ^ x \7 p \7 x + R^\7 p \7,+ 
+ V„nV„)fc, 



v Q v A v^v A h 



f \7 a a\7 f3 )h al3 . (A3) 



Using these relations, we can rewrite the bilinear expan- 
sions in a more useful form 



(V=9l% va fi) = V=gh» v {Sp^ aP a 2 +g vfj R a P x ^2V p V x -4V x V p ) - R ufjXp (V x V a - 2V a V x ) + 

+ 3R pa „ fJ a + 5R pa „ fJ \7p\7 p - g„f 3 R pa {V P V p + V P V M ) - g vf) R pp (\7 a \7 p + VV a )- 

- 2g uf3 V p DV a - 2g vfj R pa n + V^7 v V a Vp + Zg v pR^ Tp xR a TpX + , a(} R pX V p V x - 

- g^{2R pa x[iV x V p + R P xr a R pXT p) + \R 2 p<y xAg^g af 3 ~25 pv , afi )+ 

+ 2R piJ {V v V [j + V fi V v )}h ati ; (A4) 



+ ^R^VaV, - ^g vf3 R PP ,{V p V a + 3V Q V) + ^g a pR PP N p V v - R^u + 

1 1 

+ j(25 M1/ , a/ 3 - g p „g a f))R pX \7 P V \ + Rx pva Rp X + R P , P XuR a pX p - g (2S^ V t a p - g P vgaf3)R 2 p x + 

+ R pa Rufj + 2gvf3R PP R p a - g a f3R PP R p v} h al3 ; (A5) 
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(2) 



.9^/3-RV Q V M - 2it , AI „V Q ,V /3 + 2g afj R pv U + g^RVaVp + -^{S^.a/s - g^g a p)R^ - 
g^vRRafi + -^{g^uga/s — 2<5 M „ , a p)R 2 + 2Rg v pR lxa + R^R a p^ h af3 . 



(A6) 



I 

APPENDIX B: PARTICULAR RESULTS OF THE find the contribution of the commutator (38) 
CALCULATIONS IN THE BACKGROUND FIELD 
METHOD 



1 



n + 2 



- tr K^K^ — R* ua p ■ 



(Bl) 



In this Appendix one can find the results for the par- 
ticular elements of the expression (37). One can easily After a tedious algebra, we arrive at the following result: 

I 



where 



a 
b 



trU = aR 2 +bR 2 +cR 2 , 



^ ( 5xn 2 + 26xn — 2Ax — xn 3 + 6yn 

hyn 2 + lOyn - 2Ay - yn 3 + 2Axn + 8zn 



c J 2n ( y + 4x ) \ 2yn + Axn + 6zn - 2Az + hzn 2 - zn 3 
Other relevant traces are the following: 

tr (RV p P ^j = (n + 2) [ (n - l)(ai + na 2 ) -a 3 }R 2 , 



where 



also 



where 



«i 

CL2 | = 

03 



nbi +b 4 + b$-b 9 
nb 2 + b 5 + b 6 



a 4 

«5 



n(y + Ax) 



n(y + Ax) y nh _ b7 

tr [R pry V pa ^ = a 4 R\ v + a 5 R 2 , 

(n - 2)b 4 + (n 2 +n- 2)(nfe 6 + b 8 - b 9 ) + (n + 2)b 7 + 2nb 10 
(n 2 +n- 2)(6i + nb 2 + b 5 ) - (n + 2)b 3 + nb 4 - 2b w 



(B2) 



(B3) 



(B4) 



(B5) 



tr 



(v>,f - 



(3nx + ny + Ax) 2 
(y + Ax) 



16 2 { n 4 x 2 + 8 n 3 x 2 + 6 n 2 x 2 + A n 3 xy + 6 n 2 xy + y 2 n 2 + 

(y + Ax) n 



+ 10y 2 n-8y 2 -16nx 2 - 32x 2 + 32yx} R 2 ll/ + 



(y + A x) n : 



{ -24 n 3 xy - 48 n 2 xy + 32 nxy- 



- 8 z 2 n 2 + 12 z 2 n 3 + n 6 z 2 - 3 n 5 z 2 + n 4 y 2 - 11 n 3 y 2 - 128 yx - 48 n 2 x 2 + 32 y 2 - A n 4 x 2 - 36 n 3 x 2 - 



2 y 2 n 2 - 16 y 2 n - 14 n 4 zy + 128 x 2 - 2 z 2 n 4 -An 5 zx-8 n 4 zx - 32 zny + 56 zn 2 y+ 



+ 2 n 5 zy + 64 znx + 32 zn 2 x - A n 4 xy } R 2 



and finally 



tr 



ci R 2 



pi/ a (3 



c 2 Rl v + c 3 R 2 



(B6) 
(B7) 
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where the constants c\ , c 2 and C3 are given by 

12 (24 nx 2 + 8nxy + 12 n 2 x 2 - 32 x 2 + 6 n 2 a;y + y 2 n 2 ) 

ci = 2 , (138) 

(y + 4x) n 

2 

c 2 = 2 — { 24 nV + 96n 3 x 2 + 72 n 2 x 2 - 192 nx 2 + 48y 2 n- 10 yV - 128 yx- 192n:ry+ 

(y + 4x) n 2 

+ 40 n 2 xy - 32 zn 3 x + 128 znx + 32 zn 2 x - 64 zny + 16 zn 3 y + 16 zn 2 y + 52 n 3 xy + 12 n 4 xy+ 

+ 128 x 2 + 3 n 4 y 2 + n 3 y 2 + 8 z 2 n 4 + 8 z 2 n 3 + 32 y 2 } (B9) 



and 



c 3 = 2 { -4n 3 i!/ + 2n 2 xy - Anxy + 8zn 3 x - zn 3 y - Az 2 n 2 - 9z 2 n 3 + n 5 z 2 + 

(y + Ax) n 2 

+ n 3 y 2 + 32 x 2 + 16 yx - 16 n 2 x 2 - 16 y 2 - 8 n 3 x 2 - 8 y 2 n - 4 nx 2 + 3 n 4 zy -4 zny- 48 znx+ 

+ 24 zn 2 x - 24 zn 2 y + 4nz 2 + 48 zy - 96 zi } . (BIO) 
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